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Abstract 

Recently in ( Phys. Rev. D 76, 087502 (2007) and Phys. Rev. D 77, 089903(E) (2008)) a thin- 
shell wormhole has been introduced in 5-dimensional Einstein-Maxwell-Gauss-Bonnet (EMGB) 
gravity which was supported by normal matter. We wish to consider this solution and investigate 
its stability. Our analysis shows that for the Gauss-Bonnet (GB) parameter a < 0, stability regions 
form for a narrow band of finely-tuned mass and charge. For the case a > 0, we iterate once more 
that no stable, normal matter thin-shell wormhole exists. 
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I. INTRODUCTION 



Whenever the agenda is about wormholes exotic matter (i.e. matter violating the energy 
conditions) continues to occupy a major issue in general relativity It is a fact that 
Einstein's equations admit wormhole solutions that require such matter for its maintenance. 
In quantum theory temporary violation of energy conditions is permissible but in classical 
physics this can hardly be justified. One way to minimize such exotic matter, even if we 
can not ignore it completely, is to concentrate it on a thin-shell. This seemed feasible, 
because general relativity admits such thin-shell solutions and by employing these shells 
at the throat region may provide the necessary repulsion to support the wormhole against 
collapse. The ultimate aim of course, is to get rid of exotic matter completely, no matter 
how small. In the 4-dimensional (4D) general relativity with a cosmological term, however, 
such a dream never turned into reality. For this reason the next search should naturally 
cover extensions of general relativity to higher dimensions and with additional structures. 
One such possibility that received a great deal of attention in recent years, for a number of 
reasons, is the Gauss-Bonnet (GB) extension of general relativity |2|. In the brane- world 
scenario our universe is modelled as a brane in a 5D bulk universe in which the higher order 
curvature terms, and therefore the GB gravity comes in naturally. Einstein- Gauss-Bonnet 
(EGB) gravity, with additional sources such as Maxwell, Yang-Mills, dilaton etc. has already 
been investigated extensively in the literature [Sj. Not to mention, all these theories admit 
black hole, wormhole jj] and other physically interesting solutions. As it is the usual trend 
in theoretical physics, each new parameter invokes new hopes and from that token, the GB 
parameter a does the same. Although the case a > 0, has been exalted much more than 
the case a < in EGB gravity so far sl (and references cited therein), it turns out here in 
the stable, normal matter thin-shell wormholes that the latter comes first time to the fore. 

Construction and maintenance of thin-shell wormholes has been the subject of a large 
literature, so that we shall provide only a brief review here. Instead, one class jol that 
made use of non-exotic matter for its maintenance attracted our interest and we intend to 
analyze its stability in this paper. This is the 5D thin-shell solution of Einstein-Maxwell- 
Gauss-Bonnet (EMGB) gravity, whose radius is identified with the minimum radius of the 
wormhole. For this purpose we employ radial, linear perturbations to cast the motion into a 
potential-well problem in the background. In doing this, a reasonable assumption employed 
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routinely, which is adopted here also, is to relate pressure and energy density by a linear 
expression Q]. For special choices of parameters we obtain islands of stability for such 
wormholes. To this end, we make use of numerical computation and plotting since the 
problem involves highly intricate functions for an analytical treatment. 

The paper is organized as follows. In Sec. II the 5D EMGB thin-shell wormhole formalism 
has been reviewed briefly. We perturb the wormhole through radial linear perturbation and 
cast the problem into a potential-well problem in Sec. III. In Sec. IV we impose constraint 
conditions on parameters to determine possible stable regions through numerical analysis. 
The paper ends with Conclusion which appears in Sec. V. 

II. A BRIEF REVIEW OF 5D EMGB THIN-SHELLS 

The action of EMGB gravity in 5D (without cosmological constant, i.e. A = 0) is 

S = ^j VWld'x (^R + aCcB - If^.uF'^''^ (1) 

in which k is related to the 5D Newton constant and a is the GB parameter. Beside the 
Maxwell Lagrangian the GB Lagrangian Cgb consists of the quadratic scalar invariants in 
the combination 

£gB = R^- ^R.uR"" + R.upaR""'"' (2) 

in which R =scalar curvature, i?^;^ =Ricci tensor and Rfj,upa =Riemann tensor. Variational 
principle of 5* with respect to (7^^ yields 

G^, + 2aH^^ = (3) 

where the Lovelock (-ff^y) and Maxwell (T^i.) tensors respectively are 



TfMu — Ffj_aF" — -g^i,Fai3F"^ . (5) 
The Einstein tensor G^jy is to be found from our metric ansatz 

ds^ = -f (r) dt^ + jj^ + [de^ + sin^ 9 [dcj? + sin^ cjydi}^) ) , (6) 
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in which / (r) will be determined from (3). A thin-shell wormhole is constructed in EMGB 
theory as follows. Two copies of the spacetime are chosen from which the regions 



Mi,2 = {ri,2 <a, a>rh} (7) 

are removed. We note that a will be identified in the sequel as the radius of the thin-shell 
and r/i stands for the event horizon radius. (Note that our notation a corresponds to b in 
Ref. 



61]. Other notations all agree with those in Ref. [61]). The boundary, time-like surface 
Si_2 of each Mi 2, accordingly will be 

Si,2 = {ri,2 = a, a> Th}. (8) 

Next, these surfaces are identified on r = a with a surface energy- momentum of a thin-shell 
such that geodesic completeness holds. Following the Darmois-Israel formalism js] in terms 
of the original coordinates x"^ = {t, r, 9, 0, ip) , we define = (r, 9, (j), ip), with r the proper 
time. The GB extension of the thin-shell EM theory requires further modifications. This 
entails the generalized Darmois-Israel boundary conditions [9], where the surface energy- 
momentum tensor is expressed by =diag(cr,pe,p^,p^). We are interested in the thin-shell 
geometry whose radius is assumed a function of r, so that the hypersurface becomes 

S:/(r,r)=r-a(r) = 0. (9) 

The generalized Darmois-Israel conditions on S take the form 

2 {Kab - Khab) + 4a (3 Jafe - Jhab + 2PacdbK''') = -n^Sab, (10) 

where a bracket implies a jump across S, and hab = gab ~ ^aUb is the induced metric on S 
with normal vector Ua- Kab is the extrinsic curvature (with trace K)^ defined by 

ab c ^Q^aQ^b^ mng^a Q^h)^^^' K ) 

The remaining expressions are as follows. The divergence-free part of the Riemann tensor 
Pabcd and the tensor Jab (with trace J) are given by 

Pabcd = Rabcd + {Rbchda — Rbdhca) ~ (Rachdb — Radhcb) + -^R {hachdb ~ hadhcb) , (12) 
Jab = I- [2KKacKl + K,,K'''Kab - 2KacK'''Kab ' K^Kab] • (13) 
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The EMGB solution that will be employed as a thin-shell solution with a normal matter 
is given by (with A = 0) 



(14) 



with constants, M =mass and Q =charge. For a black hole solution the inner (r_) and 
event horizons (r+ = Th) are 



\ 



M 

a ± 

vr 



1/2 



(15) 



By employing the solution (14) we determine the surface energy-momentum on the thin- 
shell, which will play the major role in the perturbation. We shall address this problem in 
the next section. 



III. RADIAL, LINEAR PERTURBATION OF THE THIN-SHELL WORMHOLE 
WITH NORMAL MATTER 



In order to study the radial perturbations of the wormhole we take the throat radius as 
a function of the proper time, i.e., a = a (r). Based on the generalized Birkhoff theorem, for 
r > a (r) the geometry will be given still by (6). For the metric function / (r) given in (14) 
one finds the energy density and pressures as (6| 
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--^(^'-3(l + a^)) 



Sf) 
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p = — 
^ Air 



— + — - — M - — (1 + a^) - 2aA 



(16) 
(17) 



where I = a + f [a) /2 and A = a// (a) + a? in which 

/(« 



1 + 



Aa 



1 




(18) 



We note that in our notation a 'dot' denotes derivative with respect to the proper time r 
and a 'prime' implies differentiation with respect to the argument of the function. By a 
simple substitution one can show that, the conservation equation 



(19) 



is satisfied. The static configuration of radius oq has the following density and pressures 

V/ («o) 



(^0 



An 



- - ^ (/ («o) - 3) 
2 , f (gp) 2«r(ao) , 



(20) 
(21) 



_ao 2/(ao) flo / («o) 
In what follows we shall study small radial perturbations around the radius of equilibrium 
Oq. To this end we adapt a linear relation between p and cr as [7| 



(22) 



Here since we are interested in the wormholes which are supported by normal matter, is 
the speed of sound. By virtue of Eq.s (19) and (22) we find the energy density in the form 



a [a) 



(3^ + l)\a) 13^ + 1' 



(23) 



This, together with (16) lead us to the equation of motion for the radius of the throat, which 
reads 

3(/(a)-3-2a^)Ur^^ (^^ ^ 



V/ (a) + cfl 



An 



3 Aa ,^ 
T (/(a) 



r + 1 



a / 



After some manipulation this can be cast into 



a' + V ia) 



(24) 



(25) 



where 



V{a) = f{a)- 



+ - A 



1/3 



B 



+ B'^-A] 



1/3 



in which the functions A and B are 

,3 
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(26) 



(27) 
(28) 



We notice that V (a) , and more tediously V (a) , both vanish at a = Oq. The stability 
requirement for equilibrium reduces therefore to the determination of V"{ao) > 0, and it is 
needless to add that, V (a) is complicated enough for an immediate analytical result. For 
this reason we shall proceed through numerical calculation to see whether stability regions/ 
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islands develop or not. Since the hopes for obtaining thin-shell wormholes with normal 

n 

matter when a > 0, have already been dashed [5[, we shall investigate here only the case for 
a < 0. 

In order to analyze the behavior of V (a) (and its second derivative) we introduce new 
parameterization as follows 



16M 



m 



8Q' 



-aao, po 



-apo 



Accordingly, our new variables / (a) , ao, Po, A and B take the following forms 



and 



Po 
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ao 2/ (ao) % / (ao) 



(29) 
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(32) 
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(33) 
(34) 



Following this parametrization our Eq. (25) takes the form 

(§)%.,a,^o, 



where 



V(a) 



Via) 



a 



(35) 



(36) 



In the next section we explore all possible constraints on our parameters that must satisfy 
to materialize a stable, normal matter wormhole through the requirement V" (a) > 0. 

IV. CONSTRAINTS VERSUS FINELY-TUNED PARAMETERS AND SECOND 
DERIVATIVE PLOTS OF THE POTENTIAL 



i) Starting from the metric function we must have 
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ii) In the potential, the reahty condition requires also that 

+ 5^ > 0. 



(38) 



At the location of the throat this amounts to 



Tca: 
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or after some manipulation it yields 



an 



/(ao)-2+^<0. 



This is equivalent to 
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(40) 



(41) 



in) Our last constraint condition concerns, the positivity of the energy density, which 
means that 

5-0 > 0. 



This implies, from (31) that 



or equivalently 



|- + A(/(5„)_3) 



0<1-- + ^<4U2 
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(42) 



(43) 



(44) 



(2q y (Xq 

It is remarkable to observe now that the foregoing constraints {i — Hi) on our parameters 
can all be expressed as a single constraint condition, namely 

We plot V" (a) from (26) for various fixed values of mass and charge, as a projection into 
the plane with coordinates (3 and Qq. In other words, we search and identify the regions for 
which V" (a) > 0, in 3— dimensional figures considered as a projection in the (^, Oq) plane. 
The metric function / (r) and energy density ctq > 0, behavior also are given in Fig.s 1-4. 
It is evident from Fig.s 1-4 that for increasing charge the stability regions shrink to smaller 
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domains and tends ultimately to disappear completely. For smaller oq bounds we obtain 
fluctuations in V" (a) , which is smooth otherwise. 

In each plot it is observed that the maximum of V" (a) occurs at the right-below corner 
(say, at amax) which decreases to the left (with So) and in the upward direction (with 
Beyond certain limit (say ctmin), the region of instability takes the start. The proper time 
domain of stability can be computed from (35) as 

From a distant observer's point of view the timespan At can be found by using the radial 
geodesies Lagrangian which admits the energy integral 

/ =Eo = const. (47) 

This gives the lifetime of each stability region determined by 

1 r-- dd 

At = — , (48) 

Once Omin (amax) are found numerically, assuming that no zeros of / (a) and V (a) occurs 
for amin < a < amax, the lifespan of each stability island can be determined. We must admit 
that the mathematical complexity discouraged us to search for possible metastable region 
that may be triggered by employing a semi-classical treatment. 



V. CONCLUSION 



Our numerical analysis shows that for a < 0, and speciflc ranges of mass and charge 
the 5D EMGB thin-shell wormholes with normal matter can be made stable against linear, 
radial perturbations. The fact that for a > there is no such wormholes is well-known. 
The magnitude of a is irrelevant to the stability analysis. This reflects the universality of 
wormholes in parallel with black holes, i.e., the fact that they arise at each scale. Stable 
regions develop for each set of flnely-tuned parameters which determine the lifespan of each 
such region. Beyond those regions instability takes the start. Our study concerns entirely 
the exact EMGB gravity solution given in Ref. jsl. It is our belief that beside EMGB theory 
in different theories also such stable, normal-matter wormholes are abound, which will be 
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our next venture in this line of research. 
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Figure Captions: 

Fig. 1: V" (a) > region (m = 0.5, q — 1.0) for various ranges of /3 and Sq. The lower 
and upper limits of the parameters are evident in the figure. The metric function / (f) and 
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a"o > 0, are also indicated in the smaller figures. 

Fig. 2: V" (a) > plot for m — 1.0, q — 1.5. The stability region is seen clearly to shrink 
with the increasing charge. This effect reflects also to the ctq > 0, behavior. 

Fig. 3: The stability region for m = 1.0, q = 2.0, is seen to shift outward and get smaller. 

Fig. 4: For fixed mass m = 1.0 but increased charge q = 2.5 it is clearly seen that the 
stability region and the associated energy density both get further reduced. 
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